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Thermoelectric  (TE)  devices  can  provide  clean  energy  conversion  and  are  environmentally  friendly; 
however,  little  research  has  been  published  on  the  optimal  design  of  air-cooling  systems  for  thermo¬ 
electric  generators  (TEGs).  The  present  study  investigates  the  performance  of  a  TEG  combined  with  an 
air-cooling  system  designed  using  two-stage  optimization.  An  analytical  method  is  used  to  model  the 
heat  transfer  of  the  heat  sink  and  a  numerical  method  with  a  finite  element  scheme  is  employed  to 
predict  the  performance  of  the  TEG.  In  the  first-stage  optimization,  the  optimal  fin  spacing  for  a  given 
heat  sink  geometry  is  obtained  in  accordance  with  the  analytical  method.  In  the  second-stage  optimi¬ 
zation,  called  compromise  programming,  decreasing  the  length  of  the  heat  sink  by  increasing  its  frontal 
area  is  the  recommended  design  approach.  Using  the  obtained  compromise  point,  though  the 

heat  sink  efficiency  is  reduced  by  20.93%  compared  to  that  without  the  optimal  design,  the  TEG  output 
power  density  is  increased  by  88.70%.  It  is  thus  recommended  for  the  design  of  the  heat  sink.  Moreover, 
the  TEG  power  density  can  be  further  improved  by  scaling-down  the  TEG  when  the  heat  sink  length  is 
below  14.5  mm. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Over  the  last  several  decades,  there  has  been  a  dramatic  prog¬ 
ress  in  the  development  of  green  energy  technology  which  can 
reduce  greenhouse  gas  emissions  and  fossil  fuel  usage.  Thermo¬ 
electric  (TE)  devices,  which  consist  of  p-type  and  n-type  semi¬ 
conductors,  can  be  considered  as  a  useful  tool  to  practice  the  green 
energy  technology.  TE  devices  can  be  divided  into  two  types, 
namely  thermoelectric  coolers  (TECs)  [1-4]  and  thermoelectric 
generators  (TEGs)  [5-8].  TECs  convert  electricity  into  thermal 
energy  for  cooling  via  the  Peltier  effect,  whereas,  TEGs  convert 
thermal  energy,  say,  waste  heat,  into  electrical  power  via  the  See- 
beck  effect. 

Unlike  convectional  heat  engines  or  compression  refrigerators, 
TE  devices  are  solid-state;  they  contain  neither  moving  parts  nor 
refrigerants  [9].  Therefore,  the  whole  system  can  be  simplified  and 
operated  over  an  extended  period  of  time  without  maintenance 
[10].  TE  devices  can  produce  energy  without  using  fossil  fuel  and  can 
thus  reduce  greenhouse  gas  emissions.  However,  the  energy 
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conversion  efficiency  of  TE  devices  is  lower  than  those  of  convec¬ 
tional  heat  engines  or  refrigerators  [11  ].  The  efficiency  of  TEGs  and 
the  coefficient  of  performance  (COP)  of  TECs  are  functions  of  not 
only  the  figure  of  merit  (ZT)  but  also  the  temperature  difference 
across  the  devices  [12].  ZT  is  the  performance  index  of  a  thermo¬ 
electric  material.  Its  value  is  relatively  low  (about  1.0)  for  the  best 
existing  commercial  TE  cooling  modules  whereas  that  for  conven¬ 
tional  air-conditioning  system  is  about  4.0  [13].  Consequently, 
a  strategy  for  improving  the  performance  of  TE  devices  is  needed. 

In  reviewing  past  research  concerning  thermal  design  of  TEGs, 
a  number  of  studies  have  been  reported.  For  example,  Esarte  et  al. 
[14]  employed  a  theoretical  method  to  analyze  the  influence  of  the 
design  parameters  of  heat  exchangers  on  the  power  supplied  by 
a  TEG.  The  theoretical  results  well  matched  the  experimental  values 
for  low  flow  rates  but  not  for  high  flow  rates.  Chen  et  al.  [15]  found 
that  heat  transfer  irreversibility  affected  the  performance  of  TEG 
and  thus  had  to  be  considered  in  analysis.  A  TEG  system  combined 
with  the  heat  exchangers  at  both  the  hot  and  cold  side  was 
numerically  modeled  by  Astrain  et  al.  [16].  Their  results  showed 
that  when  the  thermal  resistances  of  heat  exchangers  on  both  sides 
of  the  TEG  were  decreased  by  10%,  the  TEG  output  power  was 
increased  by  8%. 

Recently,  waste  heat  has  been  recovered  for  further  usage. 
Waste  heat  can  be  used  for  space  and  water  heating  [17,18], 
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Nomenclature 

Too 

Fluid  inlet  temperature  (K) 

tf 

Fin  thickness  (mm) 

Ac 

Cross-sectional  area  (mm2) 

V 

Total  volume  of  heat  sink  (mm3) 

A 

Surface  area  (mm2) 

W 

Width  (mm) 

Cp 

Specific  heat  at  constant  pressure  (kj  kg-1  K-1) 

X 

Geometry  parameter  of  the  heat  sink  (mm) 

DTe 

Depth  of  TE  element  (mm) 

(*,y) 

Real  point  in  the  compromise  programming 

Dg 

Fin-to-fin  spacing  (mm) 

(x*,y*) 

Ideal  point  in  the  compromise  programming 

E 

Electric  field  intensity  vector  (V  m_1) 

ZT 

Dimensionless  TE  figure  of  merit 

f 

Distance  function  in  the  compromise  programming 

G 

Ratio  of  the  cross-sectional  area  to  length  of  TE 

Greek  letters 

element  (mm) 

a 

Fluid  thermal  diffusivity  (m2  s-1) 

Hf 

Fin  height  (mm) 

(p 

Electric  scalar  potential  (V) 

h 

Average  heat  transfer  coefficient  of  the  fins  ( W  itT2  K-1 ) 

V 

Efficiency  (%) 

h 

Heat  transfer  coefficient  (W  m-2  K_1) 

P 

Fluid  viscosity  (N  s  m-2) 

I 

Electric  current  (A) 

V 

Fluid  kinematic  viscosity  (m2  s-1) 

J 

Electric  current  density  vector  (A  m-2) 

P 

Fluid  density  (kg  m-3) 

k 

Thermal  conductivity  (W  m-1  K_1) 

Pe 

Electrical  resistivity  (Q  m) 

L 

Length  (mm) 

Nf 

Number  of  fins 

Subscripts 

Nte 

Number  of  TEG  couple 

B 

Heat  sink  base 

P 

Output  power  of  TEG  (mW) 

base 

Base  case 

P" 

Output  power  density  of  TEG,  =P/AcJE  (mW  mnrr2) 

c 

Cold  side  of  TE  element 

Pf 

Fin  perimeter  (mm) 

eff 

Effective 

Ap 

Pressure  drop  across  the  heat  sink  (N  m-2) 

F 

Fluid 

Pr 

Prandtl  number,  =  v/a 

f 

Fin 

1 

Heat  generation  per  unit  volume  (W  m-3) 

HS 

Heat  sink 

Q 

Heat  flux  vector  (W  m“2) 

h 

Hot  side  of  TE  element 

a 

Heat  transfer  rate  (W) 

L 

External  load 

Qi 

Heat  transfer  rate  for  the  boundary  layer  flow  limit  (W) 

loss 

Heat  loss  from  the  side  surfaces  of  the  TE  element 

Qs 

Heat  transfer  rate  for  the  fully  developed  flow  limit 

max 

Maximum 

(W) 

n 

n-type  for  TE  element 

R 

Electric  resistance  (Q) 

opt 

Optimum 

S 

Seebeck  coefficient  (V  I<-1) 

P 

p-type  for  TE  element 

T 

Absolute  temperature  (I<) 

TE 

Thermoelectric  element 

Tw 

Surface  temperature  of  the  fins  (K) 

t 

Total  heat  sink  heat  transfer  area 

improving  energy  recovery  efficiency  and  system  efficiency  [19,20], 
and  enhancing  chemical  reactions  [21,22].  Moreover,  several 
studies  [23,24]  have  shown  the  promising  potential  of  using  TEGs 
for  waste  heat  recovery.  Meng  et  al.  [25]  proposed  a  TEG  model 
with  multi-irreversibilities;  they  suggested  that  the  results  could  be 
regarded  as  the  feasibility  reference  using  the  waste  heat  for  power 
generation.  Because  there  is  almost  no  cost  for  obtaining  waste 
heat,  the  low  efficiency  problem  of  TE  devices  is  not  a  critical 
issue  [24]. 

Some  studies  have  optimized  the  geometric  design  of  TE 
devices.  A  numerical  optimization  of  a  TEC  was  presented  by  Xuan 
[26].  The  results  indicated  that  the  construction  cost  of  a  TEC  was 
closely  related  to  the  cooling  power  density,  whereas  the  running 
cost  was  inversely  proportional  to  COP.  Kubo  et  al.  [27]  altered  the 
size  of  incisions  along  the  lateral  faces  of  a  TE  device;  they  found 
that  the  relationship  between  the  TE  performance  and  the  incision 
size  depended  on  the  cold  side  temperature.  Yilbas  and  Sahin  [28] 
introduced  two  parameters,  the  slenderness  ratio  and  the  external 
load  parameter,  to  analyze  the  TEG  efficiency;  their  results  showed 
that  the  higher  efficiency  could  be  obtained  for  almost  all  the 
external  load  parameters  considered  when  the  slenderness  ratio 
was  less  than  1.  Jang  et  al.  [29]  optimized  the  design  of  micro-TEGs 
using  finite  element  analysis.  High  efficiency  was  obtained  when 
the  length  of  the  thermoelements  was  large.  In  addition,  the  power 
generated  declined  with  the  cross-sectional  area  of  the  thermoel¬ 
ements,  whereas  efficiency  showed  the  opposite  trend. 


A  review  of  the  literature  shows  that  the  design  of  the  geom¬ 
etry  plays  an  important  role  in  optimizing  the  performance  of 
TEGs.  However,  few  studies  have  reported  on  the  optimization  of 
geometry  design  of  TEGs  incorporated  with  air-cooling  system.  An 
air-cooling  system  combined  with  a  heat  sink  is  commonly  used 
for  dissipating  the  heat  produced  by  electronic  devices  due  to  its 
low  unit  price,  low  weight,  and  high  reliability  [30].  Accordingly, 
the  objective  of  the  present  study  is  to  investigate  the  character¬ 
istics  of  TEGs  with  an  air-cooling  design  where  a  finite  element 
method  is  used  to  predict  the  performance  of  the  TEGs.  The  effects 
of  the  heat  sink  geometry  and  TEG  dimensions  on  performance  are 
taken  into  account.  To  improve  the  performance  of  the  TEGs,  two- 
stage  optimization  is  carried  out.  Specifically,  an  analytical  method 
is  used  to  model  the  air-cooling  system,  followed  by  employing 
a  method  of  compromise  programming  to  seek  the  optimal 
performance  of  TEGs. 

2.  Mathematical  formulation  and  modeling 

Numerical  simulations  are  adopted  to  predict  the  performance 
of  TEGs.  The  physical  and  numerical  models  are  described  below. 

2  A.  Assumptions 

To  simplify  the  TEG  system,  the  following  assumptions  are 
made: 
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Front  view  Ds  tf 


Side  view 


Fig.  1.  Schematics  of  (a)  simplified  system,  (b)  thermoelectric  couple  geometry,  and  (c)  heat  sink  geometry. 


(1)  the  system  is  steady-state; 

(2)  the  existence  of  electrodes  is  ignored  and  the  contact  resis¬ 
tances  are  neglected  because  the  length  of  the  thermoelectric 
element  is  larger  than  200  pm  [31]; 


(3)  the  material  properties  of  the  thermoelectric  elements  are 
independent  of  temperature; 

(4)  the  configurations  of  the  p-type  and  n-type  elements  are 
equivalent; 
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(5)  the  TEG  module  is  constructed  from  thermoelectric  couples 
which  are  identical  and  joined  in  series; 

(6)  the  heat  sink  is  not  included  in  the  computational  domain  but 
transformed  into  a  boundary  condition  using  the  analytical 
method;  and 

(7)  the  airflow  around  the  heat  sink  is  assumed  to  be  laminar 
because  of  small  fin  space  and  low  flow  rate  typically 
encountered  in  electrical  device  cooling  [32]. 

Based  on  the  above  assumptions,  the  simplification  of  the  TEG 
module  is  illustrated  in  Fig.  la.  One  thermoelectric  couple 
comprising  a  p-type  and  an  n-type  elements  is  simulated  in  this 
study.  In  Fig.  lb,  the  geometry  of  the  thermoelectric  element  is 
specified  by  the  element  depth  DTE,  element  width  WTE,  and  element 
length  Lte-  In  the  simulation,  the  p-type  and  n-type  elements  are 
connected  electrically  in  series  and  thermally  in  parallel. 

2.2.  Governing  equations  for  TEG 

In  order  to  analyze  the  thermoelectric  system,  the  governing 
equations  include  the  thermal,  electrical,  and  thermoelectric 
effects.  Therefore,  the  conservation  principles  of  energy  and  current 
are  considered  simultaneously  [29]: 


employed  as  the  coolant.  The  direction  of  the  airflow  is  assumed  to 
be  parallel  to  the  heat  sink  base.  Due  to  the  small  fin  spacing  and 
low  airflow  rates,  the  airflow  in  typical  plate  fin  heat  sinks  applied 
to  electronic  modules  can  be  regarded  as  a  laminar  flow  [32]. 
Therefore,  the  laminar-based  analytical  method  proposed  by  Mereu 
et  al.  [35]  is  used  to  analyze  and  optimize  the  heat  sink  system  [36]. 

In  this  method,  the  scale  of  the  total  heat  transfer  rate  dissipated 
by  the  fins  is  divided  into  two  extreme  limits  with  one  for  narrow 
channels  (Dg/LH s^O)  and  the  other  for  large  channels  (Dg/LH s-> 00 ) 
[34].  For  the  limit  of  narrow  channels,  the  airflow  is  considered  as 
a  fully  developed  flow.  The  total  heat  transfer  rate  can  thus  be 
obtained  using  [35]: 


Qs  = 


pWHSH/  Df  Ap 
1  +tf/Dgl2fiLHS Cp(  w 


Too) 


(8) 


where  p  is  the  fluid  density,  g  is  the  fluid  viscosity,  Cp  is  the  specific 
heat  at  constant  pressure,  Ap  is  the  pressure  drop  across  the  heat 
sink,  Tw  is  the  maximum  surface  temperature  of  the  fins,  and  Too  is 
the  fluid  inlet  temperature.  It  should  be  noted  that  the  fin  surfaces 
are  assumed  to  be  isothermal  (i.e.,  Tw)  for  simplicity.  For  the  limit  of 
large  channels,  the  fin-to-fin  spacing  is  large  enough  so  that  the 
airflow  is  approximated  as  a  boundary  layer  flow.  Consequently,  the 
total  heat  transfer  rate  is  expressed  as  [35]: 


V-q  =  q  (1) 

V-J  =0  (2) 

where  q,  J  ,  and  q  represent  the  heat  flux  vector,  current  density 
vector,  and  heat  generation,  respectively.  The  above  equations  can 
be  coupled  by  the  following  constitutive  equations: 

q  =  SteTJ  - /<teVT  (3) 


a  -  '  (9) 

where  kf  is  the  fluid  thermal  conductivity,  v  is  the  fluid  kinematic 
viscosity,  and  Pr  is  the  Prandtl  number. 

For  a  given  heat  sink  volume  Whs  x  Ths  x  H/,  there  is  an  optimal 
fin-to-fin  spacing  Dg<0 pt  corresponding  to  the  maximum  heat 
transfer  rate.  The  value  of  Dg>op t  can  be  obtained  by  intersecting  the 
two  asymptotes  from  Eq.  (8)  and  Eq.  (9),  which  yields: 


J  =-L(£-SteVT)  (4) 

PeJE 

where  /<tE,  StE,  and  pe, te  are  the  thermal  conductivity,  Seebeck 
coefficient,  and  electrical  resistivity  of  a  thermoelectric  element, 
respectively.  Moreover,  E  is  the  electric  field  intensity  vector  which 
is  derived  from  an  electric  scalar  potential  (p: 

E  =  -Vcp  (5) 

By  substituting  Eqs.  (3)— (5)  into  Eq.  (1)  and  Eq.  (2),  the  coupled 
governing  equations  for  the  electric  potential  and  temperature  are 
given  by: 

V-(5teT7)  -  V-(/<teVT)  =  q  (6) 

V-(— !— V<p  )  +V-(SteVT)  =  0  (7) 

\PeJE  J 


%,opt  =  2  73 ( 

LHs  y  p,a  J 


(10) 


Hence,  the  maximum  heat  transfer  rate  Qmax  corresponding  to 
Dgf0pt  can  be  obtained  by  substituting  Eq.  (10)  into  Eq.  (8)  or  Eq.  (9). 
Furthermore,  in  order  to  include  the  effect  of  the  heat  sink  on  the 
performance  of  the  TEG,  an  effective  heat  transfer  coefficient  is 
used  in  this  study  [23,32].  The  effective  heat  transfer  coefficient  is 
defined  as: 


heff 


Qmax 

Ab(TW  —  Too) 


(11) 


where  AB  is  the  base  area  of  the  heat  sink.  If  the  thermal  resistance 
of  the  heat  sink  base  is  negligible,  the  influence  of  the  heat  sink 
geometry  or  fluid  operating  conditions  can  be  characterized  by  the 
effective  heat  transfer  coefficient,  which  can  be  imposed  as 
a  boundary  condition  in  the  simulation  of  the  TEG. 


where  the  heat  generation  q  in  Eq.  (6)  generally  includes  the  2.4.  Boundary  conditions 
electric  power  J  •  E  spent  on  Joule  heating  and  on  work  against  the 

Seebeck  field  StfVT  [331.  The  boundary  conditions  for  the  present  simulation  are  as 

follows. 


2.3.  Analytical  modeling  of  heat  sink 

The  heat  sink  geometry  is  shown  in  Fig.  lc.  It  is  the  most 
common  configuration  used  in  current  applications  [34].  In  this 
figure,  the  heat  sink  comprises  a  series  of  parallel  fins  with  height 
Hf,  length  LHS,  and  thickness  tf.  Each  fin  is  spaced  by  a  gap  Dg  and 
mounted  on  the  heat  sink  base  with  an  area  of  Lhs  x  Whs-  Air  is 


Table  1 

Properties  of  thermoelectric  elements  [8]. 


Material 

Electric  resistivity 

Thermal  conductivity 

Seebeck  coefficient 

(Qm) 

(W/mK) 

(V/K) 

p-type 

1.447  x  1CT5 

1.52 

226.8  x  1CT6 

n-type 

1.447  x  lO"5 

1.52 

-226.8  x  1CT6 
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Table  2 

Comparisons  of  numerical  and  analytical  results  (Th  =  423  I<  and  Tc  =  303  K). 


Quantity 

Analytical 
results  [8] 

Grid  system  1 
(2432  grids) 

Grid  system  2 
(4600  grids) 

Grid  system  3 
(12,544  grids) 

Qh  (W) 

81.3 

81.31 

81.31 

81.31 

P(W) 

3.98 

3.966 

3.966 

3.966 

/(A) 

1.08 

1.079 

1.079 

1.079 

VTE  (%) 

4.89 

4.878 

4.878 

4.878 

(1)  For  the  numerical  validations,  constant  temperature  Th  and  Tc 
are  applied  at  the  hot  side  and  cold  side  of  the  TEG  elements, 
respectively.  Furthermore,  heat  loss  along  the  side  surfaces  of 
the  TEG  elements  is  taken  into  account  by  means  of  a  uniform 
heat  transfer  coefficient  hioss- 

(2)  When  the  heat  sink  is  considered,  an  effective  heat  transfer 
coefficient  heff  instead  of  Tc  is  employed  at  the  cold  side  of  the 
TEG  elements.  The  heat  loss  is  ignored  in  this  case. 

(3)  The  ground  voltage  is  given  as  zero  at  the  cold  side  of  the  p- 
type  element. 


2.5.  Numerical  method  and  validation 

The  TEG  model  was  built  using  ANSYS  12.0.1  software.  The 
discretization  of  the  governing  equations  and  the  finite  element 
formulation  were  obtained  using  the  Galerkin  method  [37].  When 
the  simulation  was  finished,  the  results  of  the  TEG  performance, 
including  the  heat  supplied  to  the  hot  side  Qh,  the  generated 
current  /,  the  output  power  P,  and  the  efficiency  t^te  =  P/Qh,  were 
obtained. 

The  validation  of  the  proposed  TEG  model  was  performed,  and 
the  results  were  compared  with  those  of  Chen  et  al.  [8].  The  TE 
geometry  was  chosen  to  be  consistent  with  Chen’s;  that  is,  the 
values  of  1.6  mm,  1.4  mm,  and  1.4  mm  corresponding  to  the  TE 
element  length,  width,  and  depth,  respectively.  The  external  load 
resistance  Rl  was  set  to  0.0268  Q.  The  assumptions  of  Sp  =  -Sn, 
kp  =  kn,  and  peiP  =  pe<n ,  which  are  considered  reasonable  [25],  were 
made  to  reduce  the  number  of  variables.  The  obtained  material 
properties  are  listed  in  Table  1. 

The  simulation  was  conducted  using  the  above  geometry  and 
material  properties  as  well  as  the  boundary  conditions  mentioned 


Table  3 

Properties  of  heat  sink  system  [40]. 

Material 

Property 

Aluminum  (for  fin) 

kf=  237  W/m  K 

Air  (for  coolant) 

kg  =  26.3  x  10  3  W/m  K 
p  =  1.1614  kg/m3 
p  =  184.6  x  10"7  N  s/m2 
a  =  22.5  x  10-6  m2/s 

Cp  =  1007  J/kg  I< 
v  =  15.89  x  10-6  m2/s 

earlier.  The  case  without  heat  loss  was  first  examined.  Orthogonal 
grids  were  used  for  the  simulation.  Three  grid  systems  of  2,432, 
4,600,  and  12,544  grids  were  individually  tested  to  check  the  grid 
independence.  According  to  assumption  (5),  heat  supplied  to  the  hot 
side  Qh  and  the  output  power  P  of  the  TEG  module  can  be  calculated 
from  the  corresponding  results  of  the  TEG  model  by  multiplying 
them  by  the  number  of  thermoelectric  couples.  The  numerical 
results  and  the  analytical  results  are  summarized  in  Table  2.  Good 
agreement  can  be  found  between  the  numerical  and  analytical 
results.  In  addition,  the  obtained  physical  values  in  terms  of  the 
three  grid  systems  used  are  equivalent,  as  shown  in  Table  2. 
Therefore,  the  second  grid  system  (4600)  is  adopted  for  simulations. 

The  case  with  heat  loss  was  examined  next;  the  results  are 
shown  in  Fig.  2.  In  the  figure,  Qh  and  tjTe  obtained  in  the  present 
study  for  various  heat  transfer  coefficients  hioss  are  compared  with 
those  reported  by  Chen  et  al.  [8].  The  trends  of  Qh  and  tjte  are  close 
to  those  of  Chen’s.  The  maximum  relative  errors  between  the 
numerical  results  and  the  literature  data  are  less  than  10%, 
revealing  that  validity  of  the  present  prediction.  The  slight  devia¬ 
tion  may  be  caused  by  the  difference  between  the  derived  and 
actual  material  properties,  which  are  temperature-dependent. 


3.  Results  and  discussion 

3.1.  Influence  of  heat  sink  geometry  (the  first-stage  optimization) 

In  order  to  investigate  the  influence  of  heat  sink  geometry,  the 
heat  sink  configuration  from  Loh  et  al.  [38]  is  selected.  Because  their 
heat  sink  was  used  to  cool  a  thermoelectric  module,  it  is  suitable  for 


Fig.  2.  Comparisons  of  TEG  performance  between  the  numerical  results  and  literature 

data  ( Th  =  423  K  and  Tc  =  303  I<).  Fig.  3.  Development  of  the  effective  heat  transfer  coefficient  for  various  values  of  Dg. 
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this  study.  The  base  case  of  the  geometry  had  heat  sink  length 
Lhs  =  20  mm,  heat  sink  width  Whs  =  40  mm,  fin  height  Hf=  14  mm, 
and  fin  thickness  t/=  0.2  mm.  In  addition,  the  pressure  drop  across 
the  heat  sink  was  fixed  at  40  Pa  [39].  The  material  properties 
applied  to  the  heat  sink  are  listed  in  Table  3  [40].  The  TEG  geometry 
was  1  mm  (Wje)  x  1  mm  (Dje)  x  1  mm  (Lie).  Consequently,  the 
development  of  the  effective  heat  transfer  coefficient  of  the  heat 
sink  for  various  fin  spacings  can  be  obtained  from  the  analytical 
method  mentioned  earlier;  the  results  are  shown  in  Fig.  3.  In  the 
figure,  the  two  curves  corresponding  to  the  two  extreme  limits 
intersect  at  a  point  which  identifies  the  optimal  fin  spacing  Dg>0 pt 
and  the  corresponding  effective  heat  transfer  coefficient.  Further 
comparisons  of  the  TEG  performance,  namely  the  output  power 
and  efficiency,  can  be  found  in  Fig.  4a  and  b,  respectively,  where  the 
cases  for  0.5DgiOpt,  Dgi0pt,  and  2Dgopt  are  shown.  The  effective  heat 
transfer  coefficient  of  0.5Dgopt  and  2 Dgt0pt  can  be  derived  from  Eq. 
(8)  and  Eq.  (9),  respectively.  As  shown  in  Fig.  4,  there  is  a  maximum 


value  in  each  curve  (the  electrical  current  was  varied  by  altering  the 
external  loading).  In  addition,  the  case  of  Dgi0pt  covers  the  widest 
range  of  the  generated  current  and  offers  the  maximum  output 
power  and  efficiency,  as  illustrated  in  Fig.  4a  and  b,  respectively. 
Therefore,  the  results  elucidate  the  importance  of  the  heat  sink 
optimization  on  TEG  performance,  as  reported  by  Astrain  et  al.  [16]. 
It  is  noteworthy  that  the  optimal  values  of  the  current  corre¬ 
sponding  to  the  maximum  output  power  and  the  maximum  effi¬ 
ciency  do  not  coincide. 

In  order  to  investigate  the  influence  of  external  loading,  the 
output  powers  at  various  resistance  ratios  of  external  loading  to  the 
TEG  couple  are  shown  in  Fig.  5.  Instead  of  a  parabolic  curve,  the 
curves  are  more  like  Gaussian  distributions  in  a  log-log  plot.  The 
optimal  resistance  ratios  corresponding  to  the  maximum  output 
power  deviate  from  one  slightly;  that  is,  if  the  thermal  resistance  of 
the  heat  sink  is  considered  at  the  cold  side  of  the  TEG  couple,  an 
optimal  external  resistance  larger  than  the  TEG  internal  resistance 
is  needed  to  achieve  the  maximum  output  power.  The  optimal 
external  resistance  is  equal  to  the  TEG  internal  resistance  only 
when  the  thermal  resistance  of  the  heat  sink  reduces  to  zero.  For 
simplicity,  only  the  maximum  power  or  power  density  is  consid¬ 
ered  in  the  following  discussion. 

After  the  optimal  fin  spacing  of  the  heat  sink  was  found,  the 
other  geometry  parameters  were  varied.  The  procedure  was  as 
follows:  (1)  set  the  heat  sink  geometry;  (2)  determine  the  optimal 
fin  spacing;  (3)  calculate  the  effective  heat  transfer  coefficient  of 
the  heat  sink;  and  (4)  evaluate  the  performance  of  the  TEG  and  the 
heat  sink.  In  the  following  section,  the  output  power  density  is 
discussed.  Moreover,  the  heat  sink  efficiency  is  considered  as  the 
heat  sink  performance;  it  can  be  determined  as  [40]: 

Vns  =  1~  Nf^f  0  -  If)  (12) 

where  Nf  is  the  number  of  fins,  Af  is  the  surface  area  of  each  fin,  At  is 
the  total  heat  sink  heat  transfer  area,  and  rjf  is  the  fin  efficiency, 
which  can  be  obtained  as: 

tan  hfmHf 

=  mHf 


Fig.  5.  Effect  of  external  loading  on  TEG  output  power. 
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Table  4 

Dimensions  of  heat  sink  for  each  case. 


Parameters 

V  (mm3) 

Lhs  (mm) 

Whs  (mm) 

H/(mm) 

t/(mm) 

Base  case 

11,200 

20 

40 

14 

0.2 

Case  1-1 

10,  15,  20,  30,  40 

40V2,  40^4/3,  40,  40/^273,  40/^2 

14a/2,  14^/473,  14,  14/^273,  14/V2 

1-2 

80,160/3,  40,  80/3,  20 

14 

1-3 

40 

28,  56/3,  14,  28/3,  7 

Case  2 

20 

20,  30,  40,  60,  80 

28,  56/3,  14,  28/3,  7 

Case  3 

20 

40 

14 

0.1,0.15,  0.2,  0.3,  0.4 

where  m  is  defined  as: 


777 


hPf 


\  kfAc,  J 


(14) 


In  the  preceding  equation,  Pf  represents  the  fin  perimeter,  ACf  is 
the  fin  cross-sectional  area,  kf  is  the  fin  thermal  conductivity,  and  h 
is  the  average  heat  transfer  coefficient  obtained  from  the  analytical 
model  proposed  by  Teertstra  et  al.  [34]. 


It  is  worth  noting  that  the  overall  heat  sink  volume 
(V=WH s  x  Lhs  x  Hf)  is  fixed  and  the  geometry  parameter  X  (i.e.,  Lhs, 
WHs,  Hf,  and  tfi  is  considered  in  the  range  of  0.5Xbase<X  <  2XbaSe, 
where  Xbase  is  the  geometry  parameter  of  the  base  case.  In  this 
study,  three  cases  are  considered:  varying  Lhs  (Case  1 ),  determining 
WHs  by  varying  Hf  with  a  fixed  LHs  (Case  2),  and  varying  t/(Case  3). 
Because  the  overall  heat  sink  volume  is  fixed,  Case  1  can  be  further 
separated  into  three  sub-cases,  namely,  determining  LHs  by  varying 
WHsL//(Case  1-1),  determining  IHs  by  varying  WHs  (Case  1-2),  and 
determining  Lhs  by  varying  H/(C ase  1-3).  The  detailed  dimensions 
of  the  heat  sink  for  each  case  are  summarized  in  Table  4.  The  TEG 
and  heat  sink  performance  results  for  each  case  are  shown  in 
Figs.  6-8,  respectively.  Copeland  [41]  showed  that  shorter  heat 
sinks  had  better  performance  in  terms  of  thermal  conductance. 
Fig.  6  shows  that  the  three  sub-cases  of  Case  1  enhance  the  TEG 
output  power  density  but  decrease  the  heat  sink  efficiency  when 
Lhs  is  decreased.  Similar  trends  can  be  seen  in  Figs.  7  and  8.  In 
addition,  the  TEG  performance  range  is  large  for  Case  1,  whereas 
that  for  Case  3  is  small.  The  results  imply  that  the  effect  of  the  fin 
thickness  on  TEG  performance  is  relatively  slight.  Although  Case  1 
can  offer  larger  TEG  power  density,  the  drop  in  heat  sink  efficiency 
is  significant.  It  is  thus  difficult  to  evaluate  the  best  choice  when  the 
TEG  and  heat  sink  performance  are  considered  simultaneously. 
Therefore,  an  alternative  approach  named  compromise  program¬ 
ming  is  carried  out  in  the  following  section. 


3.2.  Compromise  programming  (second-stage  of  optimization ) 

To  determine  a  compromise  between  the  TEG  power  density 
and  heat  sink  efficiency,  a  multi-objective  optimization  method 
called  compromise  programming  is  utilized.  This  method  is 


Fig.  6.  Numerical  results  for  Case  1.  (a)  TEG  power  density  and  heat  sink  efficiency  and 
(b)  TEG  efficiency. 


Fig.  7.  TEG  and  heat  sink  performances  for  Case  2. 
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Fig.  8.  TEG  and  heat  sink  performances  for  Case  3. 


Base  case  Case  1-1  Case  1-2  Case  1-3 


a  well-known  approach  in  operation  research  and  management 
science  [42].  The  following  simple  form  is  applied  in  this  study: 


(15) 


where  (x,y)  represents  the  real  point,  (x*,y*)  is  the  ideal  point,  and  / 
is  the  distance  function.  In  practical  applications,  the  heat  sink 
efficiency  and  the  TEG  power  density  are  considered  as  x  and  y, 
respectively.  For  example,  the  compromise  programming  operation 
for  Case  1  is  shown  in  Fig.  9.  If  a  heat  sink  efficiency  of  100%  and  the 
TEG  power  density  at  the  conditions  of  Th  =  423  K  and  Tc  =  303  K 
(25.59  mW  mm-2)  are  regarded  as  the  ideal  point  (x*,y*),  the 
distance  to  each  real  point  can  be  determined  using  Eq.  (15).  The 
point  that  has  the  minimum  distance  for  Case  1-1,  called  the 
compromise  point,  is  shown  in  Fig.  10.  In  this  case,  though  the  heat 
sink  efficiency  is  reduced  by  a  factor  of  20.93%,  the  power  density 


Fig.  9.  Curves  obtained  using  compromise  programming. 


b 


Fig.  10.  TEG  performance  at  compromise  points  for  the  sub-cases  of  Case  1  obtained 
using  the  compromise  programming,  (a)  TEG  power  density  and  (b)  heat  sink 
efficiency. 


can  be  enhanced  by  88.70%  compared  to  those  of  the  base  case.  If 
this  case  is  compared  with  the  results  of  the  other  geometry 
parameters,  say,  Case  2  and  Case  3,  the  results  shown  in  Fig.  11 
indicate  that  Case  1-1,  where  LH s  is  reduced  by  increasing  WHsH/, 
is  the  best  choice  when  the  TEG  and  heat  sink  performance  are 
considered  simultaneously.  Moreover,  as  mentioned  earlier,  the 
effect  of  tf  on  the  TEG  and  heat  sink  performances  is  negligible  even 
when  the  best  compromise  point  is  obtained. 

3.3.  Scaling  effect  on  TEG  performance 

As  mentioned  above,  the  case  of  determining  LH s  by  varying 
WhsH/  is  the  recommended  approach  for  achieving  the  best 
compromise  between  the  TEG  and  heat  sink  performances. 
Therefore,  the  scaling  effect  on  the  TEG  performance  is  taken  into 
account  in  this  case  for  further  investigation.  In  order  to  keep  the 
thermal  conductance  and  electric  resistance  of  the  TE  couples 
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Fig.  11.  TEG  performance  at  compromise  points  for  three  cases  obtained  using 
compromise  programming,  (a)  TEG  power  density  and  (b)  heat  sink  efficiency. 


Fig.  12.  (a)  Scaling-down  effect  on  TEG  performance  and  (b)  scaling-up  effect  on  TEG 
performance. 


4.  Conclusions 


constant,  the  scaling  approach  with  a  constant  geometry  factor 
(^c,te/^te)  is  employed  [2].  Accordingly,  three  kinds  of  the  TEG 
geometry  (Wje,  Dje,  ^te)  with  the  same  geometry  factor  are 
considered:  the  base  case  (1  mm,  1  mm,  1  mm),  the  scaling-down 
case  (0.5  mm,  0.5  mm,  0.25  mm),  and  the  scaling-up  case  (2  mm, 
2  mm,  4  mm).  The  scaling-down  effect  on  TEG  performance  is 
shown  in  Fig.  12a  and  the  scaling-up  effect  on  TEG  performance  is 
shown  in  Fig.  12b.  In  Fig.  12a,  a  larger  output  power  density  can  be 
obtained  by  scaling-down  the  TEG  size  when  the  heat  sink  length 
Lhs  is  below  about  14.5  mm.  This  implies  that  the  drop  in  cross- 
section  area  of  the  TEG  is  smaller  than  that  in  output  power  and 
thus  a  larger  power  density  is  obtained  when  Lhs  is  below  14.5  mm. 
Therefore,  the  heat  sink  design  is  a  key  point  for  the  scaling-down 
case.  However,  the  power  density  cannot  be  improved  by  scaling- 
up  the  TEG  size  under  the  considered  range  of  Lhs  as  shown  in 
Fig.  12b.  The  lower  power  density  is  due  to  the  fact  that  the  increase 
in  output  power  is  lower  than  that  in  cross-section  area  of  the  TEG 
for  the  scaling-up  case.  The  efficiency  of  the  larger  TEG  is  better  due 
to  the  larger  temperature  difference  across  the  TEG  elements. 


A  TEG  combined  with  an  air-cooling  system  has  been  investi¬ 
gated  in  this  study  through  a  method  of  two-stage  optimization  for 
designing  heat  sink  geometry.  An  analytical  method  is  employed  to 
model  the  air-cooling  system  and  the  finite  element  method  is 
utilized  to  simulate  performance  of  the  TEG.  In  the  first-stage 
optimization,  an  analytical  method  using  the  concept  of  effective 
heat  transfer  coefficient  is  adopted  to  design  the  optimal  fin-to-fin 
spacing  of  the  heat  sink.  The  results  show  that  better  TEG  perfor¬ 
mance  can  be  achieved  after  the  first-stage  of  optimization.  The 
optimal  resistance  ratio  of  external  loading  to  the  TEG  couple 
slightly  deviates  from  one.  In  the  second-stage  optimization,  the 
compromise  programming  method,  a  compromise  between  the 
TEG  performance  and  the  heat  sink  performance  is  obtained  when 
the  heat  sink  volume  is  fixed.  The  results  suggest  that  decreasing 
the  length  of  heat  sink  by  increasing  the  frontal  area  (V\4isL//)  of  the 
heat  sink  is  the  best  approach.  In  this  case,  at  the  compromise  point, 
the  TEG  output  power  density  is  improved  by  88.70%  and  the  heat 
sink  efficiency  is  reduced  by  only  20.93%  compared  to  the  base  case. 
The  influence  of  the  fin  thickness  is  negligible  even  when  the  best 
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compromise  point  is  found.  The  power  density  of  the  TEG  can  be 
further  enhanced  by  scaling-down  the  TEG  size  with  a  constant 
geometry  factor  when  the  heat  sink  length  is  below  14.5  mm. 
Alternatively,  the  power  density  cannot  be  improved  by  scaling-up 
the  TEG  size.  However,  a  better  efficiency  is  obtained  by  scaling-up 
the  TEG  size  due  to  a  larger  temperature  difference  across  the  TEG. 
The  results  have  provided  a  useful  guideline  for  designing  TEGs 
combined  with  an  air-cooling  system. 
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